A study of nuclear stopping in central symmetric nuclear collisions at intermediate 

energies 
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Nuclear stopping has been investigated in central symmetric nuclear collisions at intermediate 
energies. Firstly, it is found that the isotropy ratio, Ri so , reaches a minimum near the Fermi 
energy and saturates or slowly increases depending on the mass of the system as the beam energy 
increases. An approximate scaling based on the size of the system is found above the Fermi energy 
suggesting the increasing role of in-medium nucleon-nucleon collisions. Secondly, the charge density 
distributions in velocity space, dZ/dv\\ and dZ/dv±, reveal a strong memory of the entrance channel 
and, as such, a sizeable nuclear transparency in the intermediate energy range. Lastly, it is shown 
that the width of the transverse velocity distribution is proportional to the beam velocity. 



The study of transport phenomena in nuclear reactions 
at intermediate energies is of major importance in the 
understanding of the fundamental properties of nuclear 
matter (for an introduction seeQ). The comparison of 
the predictions of the microscopic transport models (see 
for instance H Q H) with experimental data can 
help improve our knowledge of the basic ingredients of 
such models: namely the nuclear equation of state and 
as such, the in-medium properties of the nucleon-nucleon 
interaction. In this context, it is mandatory to test the 
different models over a large systematics in system size 
and incident energy. Among the different issues that can 
be addressed in such a framework, the question of the 
thermalization of the system (in particular the damp- 
ing of the momentum distribution) in strongly dissipa- 
tive collisions is among the most debated in view of its 
connection with the search for a phase transition of the 
liquid-gas like type@,0|- 

In this paper, we investigate this problem by study- 
ing the stopping power in nuclear reactions at interme- 
diate energies. We take advantage of the high quality of 
the data collected for a large variety of systems with the 
INDRA multidetector both at GANIL (see for instance 
H El El El El) and at GSI [HE!- Symmetric sys- 



tems with total sizes between 80 and 400 mass units and 
with incident energies between 25 and 100A MeV have 
been considered. Event selection is performed as follows: 
in view of the experimental detection thresholds in the 
backward part of the detection device, we have only con- 
sidered charged products which are emitted in the for- 
ward velocity space in the centre-of-mass of the reaction. 
In this way the effect of the thresholds are minimized as 
well as the associated distortions on the global variables 
considered in the following. The selection criterium as 
far as the quality of the data is concerned was the fol- 
lowing: events were retained if the total detected charge, 
Ztot (neutrons were not detected) in the forward hemi- 
sphere of the centre-of-mass was larger than 90% of half 
the total charge of the considered system. Contrarily to 
the study of the FOPI data presented in El, no attempt 
has been made to symmetrize the data. 

In the following, we concentrate on central collisions. 
This raises the question of the selection of those events 
on criteria which do not induce auto-correlations. Since 
we are interested in quantities involving both the trans- 
verse and longitudinal directions, it is not suitable to use 
vector variables such as, for instance, the transverse en- 
ergy [HE3| Thus, a scalar variable is needed and the 
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Figure 1: Data for Xe+Sn collisions at 45 A MeV: (a) Distri- 
bution of total charged particle multiplicity u to t ■ (b) Evolu- 
tion of the mean isotropy ratio Rise, as a function of the total 
multiplicity, vtot- (c) Bi-dimensional plot showing the corre- 
lation between the isotropy ratio and the total multiplicity. 
Different intensities (in log scale) represent different numbers 
of measured events, (d) As (c), but after renormalising the 
number of events in each bin according to the evolution of the 
cross-section with the multiplicity from (a). 



most natural one is the multiplicity of detected charged 
products, v to t- In that sense, a minimum bias selection is 
used. We now consider the isotropy ratio, Ri SO , defined 
as : 
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where E± ( En ) is the cm. transverse (parallel) en- 
ergy and the sum runs over all products with the above- 
mentioned selection. We first consider the evolution of 
Riso as a function of v to t • Another possibility was to 
consider the evolution as a function of the particle mul- 
tiplicity emitted in the forward direction but this does 
not affect the conclusions. The four panels in Figure 
illustrate the selection method of the central collisions. 
Figure QJi shows the distribution of v to t while Figure 
displays the evolution of the mean value of Ri SO as a 
function of v to t ■ The mean value of Ri SO displays a S-like 
shape and reaches an asymptotic value hereafter noted 




20 30 40 50 60 70 80 90 100 
Beam Energy (A MeV) 



Figure 2: Distribution of R^o in Xe+Sn central collisions. 
The incident energy, mean value and variance of Rlto are 
indicated in each panel. The last panel shows the evolution 
of the RMS (squares and left scale) while triangles and the 
right scale correspond to v cen defined in the text. Note that 
both scales do not start from zero. 



R^g, and its corresponding mean value < >, for 

the largest values of v to t- Strong fluctuations are ob- 
served for extreme values of v to t where the statistics is 
very low (typically less than 10 events, see Figure^). In 
the following, we calculate the asymptotic value of Ri SO 
by considering the mean value of Ri SO for events with a 
multiplicity larger than v cen (these are plotted as a func- 
tion of the incident energy as triangles (right scale) in the 
last panel of Figure EJ. v cen is defined as the multiplicity 
of charged products at which Ri SO reaches its asymptotic 
value. Figures Oh and^l show the bi-dimensional corre- 
lation between v to t and Ri SO - It is worth noting from Fig. 
Gt that the largest fluctuations of Ri SO are not associated 
with the largest values of the multiplicity. This could 
suggest that the selection of central events by means of 
vtot is inappropriate. Let us discuss briefly this point. 
A first effect to consider is the influence of the statistics 
of the multiplicity distribution : this has been accounted 
for in Fig. GJ1 by dividing the contents of each bin by the 
number of events with the corresponding multiplicity ob- 
tained from Fig. QJi- The large fluctuations observed for 
intermediate values of the multiplicity have disappeared 
demonstrating that this effect is to a large extent due to 
statistics. 

Another argument in favour of the present selection of 
central collisions concerns the velocity-dependent charge 
densities [isL fl9| (defined later in this paper). Such dis- 
tributions (not shown here except for central collisions, 
see Figure do not change anymore for multiplicities 
vtot > v cen . This shows that for multiplicities larger 
than v cen , the events have the same kinematical charac- 
teristics and that they only differ by the sorting variable, 
namely the multiplicity. We call these events "central 
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Figure 3: Top: Evolution of the mean value of ii^l™ (and 
the associated uncertainty) as a function of the beam energy 
for central collisions and for various systems: (stars) Au+Au, 
(dots) Xe+Sn, (triangles) Ni+Ni and (squares) Ar+KCl. The 
solid line is the value of Ri SO corresponding to the initial nu- 
cleoli momentum distributions consisting of two sharp Fermi 
spheres separated by a relative velocity corresponding to the 
considered beam energy. Bottom: Same as above for the 
scaled isotropy ratio: < RH" > /A 1 * 13 . The solid line is a 
linear fit to all data points above 40A MeV. 

events" : they correspond to roughly 1% of the reaction 
cross-section and are the only events considered in the 
following. 

The present study has been undertaken for a variety 
of systems in the intermediate energy range. Let us first 
consider the nearly symmetric Xe+Sn system. Figure 
|2l shows examples of -R£|™ distributions between 25 and 
100A MeV incident energy. The distributions are not 
gaussian-like and the most probable values do not dras- 
tically change and remain close to 0.5-0.6 as the beam 
energy increases. The evolution of the mean value is thus 
to a large extent governed by the evolution of the width 
of the distributions. These latter are plotted in the last 
panel (squares and left scale). A minimum is observed 
near 40-50A MeV close to the Fermi energy. Also shown 
are the values of v cen (triangles) in the last panel. They 
increase from about 22 at 25A MeV up to 47 at 100A 
MeV. 

Figure shows the evolution of < > for various 

symmetric systems whose total sizes vary from about 80 
up to 400 mass units. The upper panel of the figure shows 
that, below 40^4 MeV, the isotropy ratio is nearly inde- 
pendent of the size of the system and slowly decreases. 
However, as the beam energy crosses the Fermi energy, 
one observes an increase of the isotropy ratio depending 
on the size of the system: the larger the size, the larger 
< Rig" >■ The solid line on Figure 0-top is the value 



of the isotropy ratio in the entrance channel due to the 
Fermi distribution of the nucleons inside the projectile 
and the target. Thus, the distance between the exper- 
imental points and the solid line is a "measure" of the 
influence of the dissipation on the isotropy ratio. 

The general evolution of the data is interpreted as a 
transition from the dominant influence of the mean field 
at low energy (one-body dissipation which does not de- 
pend on the size of the system) towards the dominance 
of in medium nucleon-nucleon collisions at higher energy 
(two-body dissipation which depends on the size of the 
system). 

In this case, the key quantity is the ratio between the 
nucleon mean free path and the size of the system. This 
is tentatively put in evidence in Figure (lower panel) 
where the reduced isotropy ratio R re d —< Rlfo > /A 1 / 3 
has been plotted for all considered systems (here A has 
been arbitrarily taken as half the total mass of the sys- 
tem as we are dealing with symmetric systems). For 
medium and heavy systems (Xe+Sn and Au+Au), the 
scaling is evidenced for beam energies larger than 40-50A 
MeV while for lighter systems (Ar+KCl and Ni+Ni), it 
is more and more verified as the beam energy increases. 
These features outline the role of the size of the system 
as far as the damping of the nucleon momentum distribu- 
tion is concerned. This scaling and the rather low values 
of < R^so > suggest that the mean free path of the nu- 
cleons inside the medium is rather long at such incident 
energies. It is of the same order of magnitude or even 
larger than the size of the system. One should also have 
in mind surface effects in the sense that even for cen- 
tral collisions, those nucleons that are more localized in 
the vicinity of the surface have little chance to experi- 
ence hard collisions: All in all, this means that two-body 
dissipation is not enough to drive the system towards 
thermalization. 

Up to now, we have only considered a global variable, 
Riso, built on an event by event basis. We now consider 
a variable which is averaged over all the selected events: 
namely the velocity-dependent charge density, dZ/drj^ ed 
and dZj drf£ d respectively along the beam axis and along 
one of the transverse directions. These quantities are 
built by considering all charged particles emitted in the 
forward centre of mass direction, adding the particles 
(weighted with their respective nuclear charges) which 
are localized in the same velocity bin. These distribu- 
tions are well-suited to a comparison with microscopic 
transport models, as shown for instance in 0|. The 
same distributions have been studied for FOPI data in 
[lfj . In order to compare easily the same system at dif- 
ferent beam energies, scaled velocities have been used: 
77 = v/vcm where v cm is the projectile velocity in the 
centre-of-mass frame. 

Figure0]shows the evolution of the two densities for the 
system Xe+Sn at various incident energies between 25 
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Xe+Sn at 39 AMeV 
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Xe+Sn at 65 AMeV 
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3 Xe+Sn at 100 AMeV 
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Figure 4: Top: dZ/dr]^ ed (black points) and dZ / drff d (open points) in the parallel and in one of the two transverse directions 
for central Xe+Sn collisions for various incident energies indicated in each panel. The distributions are normalized in such way 
that the total charge corresponds approximately to half the total charge of the sytem. Last panel: Evolution as a function of 
the incident energy of the width of the charge density distribution in the transverse direction, for (stars) Au + Au, (circles) Xe 
+ Sn, (triangles) Ni + Ni and (squares) Ar + KC1 systems. 



and MeV. There is a clear correlation between Ri SO 

and the behaviour of the two charge densities. Even at 
25A MeV, the two distributions are not identical, which 
would indicate the occurrence of full stopping. However, 
this is the incident energy for which the two distribu- 
tions are the closer to each other and consequently, it 
corresponds to the largest value of Ri SO . For all the con- 
sidered energies, a strong memory of the entrance chan- 
nel is evidenced: the widths of dZ/dr]^ ed are larger than 
the widths of dZjdr\ r l d . A closer look at the transverse 
distributions shows that they are quite similar whatever 
the incident energy. This is illustrated in the last panel 
of Figure Q] The reduced width of dZj drff d shows a re- 
markable constancy as a function of the incident energy 
indicating that the transverse collective motion is pro- 
portional to the beam velocity. Moreover there is no sig- 
nificant dependence on the size of the system, the same 
constancy is observed for other systems with different 
sizes : Au + Au, Ni + Ni and Ar + KC1 (see Figure |U 
last panel) . It would be interesting to study the origin of 
this effect in microscopic transport calculations. 

To conclude, we have studied the behaviour of two 
observables, the isotropy ratio Ri SO and the velocity- 
dependent charge density dZ/drf in central symmetric 
nuclear reactions at intermediate energies. The evolu- 
tion of Riso as a function of the incident energy shows a 
minimum followed by a moderate increase for most stud- 
ied systems or a plateau around the Fermi energy for the 
lightest ones. Rather low values have been measured indi- 



cating that full stopping is far from being achieved in this 
energy regime. As far as medium-mass and heavy sys- 
tems are considered, a scaling law in terms of the size of 
the system is observed suggesting the role of in-medium 
nucleon-nucleon collisions in the slow increase of the stop- 
ping power. The study of the velocity-dependent charge 
density in the parallel and transverse directions with re- 
spect to the beam shows a strong memory of the entrance 
channel in the sense that the densities in the parallel and 
transverse directions do not coincide suggesting trans- 
parency as a key feature of nuclear reactions at interme- 
diate energies. Last, the width of the transverse charge 
density using reduced velocities is invariant in the entire 
intermediate energy range. 

Our findings concerning the stopping power in cen- 
tral nuclear collisions are in full agreement with a recent 
study based on very similar observables performed by the 
FOPI Collaboration at higher incident energies up to the 
GeV region Our results as well as those of FOPI 

should be compared with the predictions of existing mi- 
croscopic transport models. As such, we believe that they 
constitute a very severe test of these models over a wide 
range of system sizes and incident energies. 
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